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We develop an analytical theory of the localization-delocalization transition for a disordered Bose 
system, focusing on a Cooper-pair insulator. We consider a chain of small superconducting granules 
coupled via Josephson links and show that the low-temperature tunnelling transport of Cooper pairs 
is mediated by a self-generated environment of dipole excitations comprised of the same particles 
as the tunnelling charge carriers in accord with the early notion by Fleishman, Licciardello, and 
Anderson [l[. We derive an analytical expression for the current- voltage characteristic and find 
that at temperatures, T, below the the charging energy of a single junction, Ec, the dc transport 
is completely locked by Coulomb blockade effect at all voltages except for a discrete set of resonant 
ones. At T > Ec the combined action of disorder and temperature unlocks the charge transport, 
since the environment excitation spectrum becomes quasi-continuous according to a Landau-Hopf- 
like scenario of turbulence, and the conductivity acquires an Arrhenius-like thermal activation 
form. The transition from the localized to delocalized behaviour occurs at T = Ec which corresponds 
to the onset of turbulence in the spectral flow of environmental excitations with Reynolds number 
TZe = [ksT/Ec) — 1. The proposed theory breaks ground for a quantitative description of dynamic 
and quantum phase transitions in a wealth of physical systems ranging from cold atoms in optical 
lattices, through disordered films and wires to granular and nanopatterned materials. 



INTRODUCTION 

In strongly disordered systems the low-temperature 
charge transport occurs via tunnelling between localized 
states, having in general different energy levels. Thus 
tunnelling is possible only if charge carriers can either 
emit the excess or absorb the deficit in energy to accom- 
modate the difference between the initial and final states 
[see inset in Fig. [1^)] . This process is referred to as the 
energy relaxation to a bosonic environment. The most 
common environment is a phonon bath; recent studies re- 
vealed, however that in Josephson-junction arrays 
( J JAs) and granular materials another relaxation mecha- 
nism naturally arises: emission/absorption of dipole exci- 
tations (charge-'anti-charge' pairs) comprised of the same 
particles that carry the current. 

This self-generated dipole environment possesses an in- 
finite number of degrees of freedom and thus serves as the 
thermostat itself, dominating over the usual relaxation 
mechanism via phonons, which becomes inefficient at low 
temperatures. Indeed, in highly disordered materials the 
standard electron-phonon interaction via a deformation 
potential should produce phonons with wavelength on the 
order of the electron wavelength A, which, in its turn, 
is on the order of the localization length. In granular 
materials or JJAs the role of the localization length is 
taken by the granule size (superconducting granules in 
the case of a JJA). Thus the characteristic phonon en- 
ergy is Eph — hqs ~ hs/X, and the corresponding phonon 
temperature is Tph ~ Eph/ks which for typical local- 



ization lengths on the order of 10 nm is about 10 K. At 
temperatures T <^ Tph relaxation via phonons becomes 
inefficient and gives way to non-phonon mechanisms. 



Since phonons arc not involved into the process of tun- 
nel charge transfer, the related heating processes in JJAs 
and granular systems, which can be referred to as elec- 
tronic insulators, can be expected significantly reduced as 
compared to conventional insulators. The suppression of 
heating when relaxation occurs via energy exchange with 
a non-phonon environment was indeed found in single- 
junction systems 



Hi. 



To proceed to a theory of bosonic insulator we con- 
sider an exemplary system, a one-dimensional (ID) JJA, 
a tunable and experimentally accessible realization of a 
generic ID strongly interacting (charged) disordered bo- 
son chain, modelling a wealth of physical systems ranging 
from helium in vycor [l2| , localized Cooper pairs near the 
superconductor-insulator transitions in thin film [ll,[ii], 
and ultracold atoms (TBHSl- same time, in the 



limit where the capacitance of a single junction well ex- 
ceeds its capacitance to the ground, the ID JJA is equiv- 
alent to a chain of coupled quantum rotors and allows 
for an analytical description. This makes it a unique 
theoretical laboratory for studying disorder-induced non- 
conventional insulators with physics arising from the in- 
terplay between strong disorder, quantum fluctuations, 
and strong Coulomb interactions. 



MODEL AND TUNNELLING TRANSPORT IN 
ONE-DIMENSIONAL JOSEPHSON JUNCTIONS 
ARRAY 

The behaviour of JJAs is controlled by two compet- 
ing parameters, Ec = ^E^^ /N, the average charging 
energy of a single superconducting junction and Ej, the 
average Josephson coupling energy between neighbour- 
ing superconducting granules. Here N is the length of 



the array and E, 



(i) 



VC^'-* the charging energies of the 



individual junctions defined by their capacitances, C^'^ 
If Ec < Ej the system is superconducting and Ec > Ej 
corresponds to an insulating state. We focus on a deep 
insulating state, Ec 3> Ej, treating Ej as a perturbation. 
At the same we let the disorder-induced dispersion in iJc'-* 
be also large as compared to Ej, such that a continuous 
conduction band due to overlapping localized wave func- 
tions does not form. We start with a qualitative picture 
of Cooper pair transport, that can be constructed in the 
context of simple system consisting of a single supercon- 
ducting granule confined between the two electrodes. 

The tunnel current is controlled by the intensity of the 
relaxation and assumes the form [19| 

/ cx cxp(-£;/M^) , (1) 

where E is some characteristic energy associated with 
the tunnelling process and W quantifies the relaxation 
rate. One can expect W to be proportional to the den- 
sity of environmental excitations. At high temperatures 
(exceeding the local gap in the excitations spectrum) the 
equipartition theorem yields W T. At low temper- 
atures, where the discreteness of the environment spec- 
trum becomes essential, the relaxation rate is strongly 
suppressed and the tunnelling current is locked. 

Although this simple line of reasoning does not apply 
straightforwardly and neither is a single quantity W char- 
acterizing the relaxation rate well defined , it gives a 
good qualitative idea of the tunnelling transport in in- 
sulators and the nature of the localization-delocalization 
transition. To avoid the problem of boundary conditions 
at the current leads we close JJA into a ring, the current 
being induced by a time changing magnetic field per- 
pendicular to the ring plane, see Fig. [1^). The current 
flowing through the series of junctions is the same along 
the ring and one can therefore calculate it across an arbi- 
trary junction (i.e. between an arbitrary pair of adjacent 
granules). The current is determined by the correlation 
function of superconducting phases along the ring, de- 
scribing the probability of the energy exchange between 
the tunnelling Cooper pair and the environment. The 
correlation function is composed of matrix elements be- 
tween the states where the number of Cooper pairs at 
a given granule differs by one, i.e. the matrix elements 
between the states \n^'^'>) and |7i(*+^)) (ti^*) is the number 
of Cooper pairs at the i-ths junction between supercon- 
ductor i and i + 1). 



3) H(t)A-(N+1)Vt 




b) 





ground state 



excited state 



FIG. 1. Model system a) A circular array of Josephson 
junctions containing + 1 junctions (superconductors are 
dark gray, tunnel barriers blue), where ^''^ and C'-'-' are the 
phase difference and the capacitance across the junctions be- 
tween grains i and i + 1, respectively. The current in the 
array is induced by a time dependent perpendicular magnetic 
field H{t). The total voltage follows from the time deriva- 
tive of the total flux <l?(t) = H{t)A [A denotes the area of 
the ring). The inset shows a sketch of a tunnelling process 
between two mesoscopic potential well. Since the energy lev- 
els in these wells are in general different (due to structural 
disorder), quantum mechanical tunnelling is only possible if 
the tunnelling particle emits (or absorbs) bosonic excitation 
which can accommodate the energy difference 5e between the 
respective energy levels, b) Illustration of Left: ground state 
and Right: exicted current states in the JJA ring 



Switching between these states means exciting a cur- 
rent state in the ring, see Fig.[T}D), i.e. removing a Cooper 
pair from the i-ths granule and adding it to the z -t- 1- 
ths granule and so on eventually encircling the system, 
2— >z — 1— >i — 2... — >i-|-2— >-j-|-l. In the insulating 
limit the current is due to tunnelling and local Coulomb 
blockade effects make the local environmental spectrum 
discrete (with a minigap of the order of the charging en- 
ergy of a single junction Ec), which means that environ- 
mental dipole excitations get localized. Localization of 
environmental excitations suppresses the relaxation and 
thus completely locks the current. However, as we show 
below the combined effect of disorder and temperature 
gives rise to Landau-Hopf-like turbulence in the spec- 
tral flow of the environment levels as soon as the ratio 
T / Ec, which plays the role of the Reynolds number in 
our system, becomes large. As a result the spectrum at 
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temperatures above the minigap becomes locally quasi- 
continuous unlocking the current. At T Ec the spec- 
tral flow of the environmental levels can be characterized 
as rather "laminar," the local spectrum retains its dis- 
crete nature and tunnelling current remains significantly 
suppressed. 

The JJA loop is described by the quantum rotor La- 
grangian: 

C = Cc + Cj, (2) 

= i E ^^'^ (x. - x^+lr + ^ E ^o,x!, (3) 



-Cj = - ^ Ey cos {xi - Xt+i) 



(4) 



where Xi is the phase of the order parameter at the su- 
perconductor i, xn+1 = Xo, C^'' and Ej^ are the capac- 
itance and the Joscphson coupling of the junction con- 
necting the i-ths and (j + l)-ths grains, respectively - 
we use the shorthand notation {i) for junctions {i, i + 1). 
Quantities related to junctions use the superscript (i) 
and subscripts i label single grain properties. The terms 
£1:7 and Cj describe the charging energy and the Joseph- 
son coupling between adjacent superconducting islands, 
respectively. The capacitance to the ground of the i- 
ths superconducting grain is Co,i. Here the convention 
h = c = kB = iis used and lengths are measured in units 
of a single junction size. We consider the deep insulating 



state where E^J^ = e^/C'-''\ E^^ , ^ c^/Cqa > E^ . The 
bias V is the voltage drop across the single junction. We 
address the most common experimental situation where 
max Co, i <C minC^*' and focus on the case where the 
charge screening length A = (C/Cq)^/^ > N, the gener- 
alization onto the opposite limit is straightforward. The 
bar denotes averaging over grains. 




FIG. 2. Matrix-Theta function Semi-logarithmic surface 
plot of the matrix i9-function which defines the current in the 
systems as function of voltage and temperature fro N = 8, 
a = 10~*i5c, and charging energy variance a = 0.056. 



Treating Cj as a perturbation [K 
percurrent in junction (i): 



one finds the su- 



= -4e[£;yT dtsin 2eV'-'h Im[/C(t)]. (5) 



Here F*^*' is the voltage drop at the junction between 
grain i and i + 1 and 



IC{t) = (exp[*0«(i)]exp[-z0«(O)])«, 



(6) 



where 4>^^^ (t) is the phase difference across the junction, 
describing the fluctuations around the mean value deter- 
mined by the external voltage. Using equations ([5]) and 
(pl) one can derive an analytical description of Cooper 
pair tunnelling transport in ID JJAs. The Fourier trans- 
form of the phase correlation function, JC{uj), 



{n},{m} 



(7) 



where \{n}) and denote the wave functions and en- 
ergy levels for configuration {n} of Cc, accordingly. The 

Gibbs weight is given by P{7i} = Z^^e with Z = 

E^'^/Ec [this expression for Z is obtained by using Pois- 
son's summation formula]. /C(w) has the physical mean- 
ing of the probability of exchanging an energy quanta 
to with the environment, which is comprised of the ex- 
citations of the system itself. This means, the system 
generates its own environment. The charging energy of a 
junction (i) is equal to e*^*)(n'^'^) = 2£'c''' [?^'■'•']^, where 
the quantum numbers n^'^ = 0, ±1, . . ., i = 0, . . . ,N 
have a physical meaning of an excess (or deficit) num- 
ber of Cooper pairs at the junction. Therefore = 
^ e(')(n(*)), where {n} = n'^^\ . . .). The wave func- 

tion describing the quantum states of the annular JJAs 
is 



N 



(8) 



j=0 



where J2i=o 4'^^^ = and '!/'n(0) = exp{2?^0}/\/27r is the 
wave function of the quantum rotor. The matrix element 
of the operator exp[— ic/)^'^ (0)] in ([6]) is nonzero for tran- 
sitions where either (a) only ti^*' changes by one or (b) 
where n*^*) is fixed but the quantum numbers of the rest of 
the junctions change by one simultaneously. Integrating 
the most singular contributions to the current explicitly 
one arrives at 



2TNEc 



■2TNEc 



1? 



3V/E,,T)I, 



(9) 
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where X 



eI^'^/e,, Z = 



E 



{m}' 



-2_E<')mJ/T^ z = V/Ec- In the partition function 



the sum is taken over all possible combinations of in- 
teger numbers {m} = (mi, m2 . . . mjv). The generalized 
Jacobi i^-functions is defined as [20|: 



{m} 



(10) 



where a = (1,1,..., ly /N and the sum is taken over all 
integer vectors m = (mi, m2, . . . , m^r). The matrix 7~ is 
parametrized via the vector e = (e*^^\ e'-^-', . . .e^^^)"^: 

2nT 



Er 



-Eli 



2{NE,y 



Hi 



The quantity a E^ introduces the finite width of the 
quantum levels. The matrix T is positively defined, since 
the matrix H has one zero eigenvalue corresponding to 
the eigenvector Jv'^'^ ~ e while other eigenvalues of H are 
of the order of unity. 

We parametrize the dispersion in charging energy as 
6E^ = {Y.{Ec^ - EcflNfl'^ = aEc. In the "clean" 
limit, a = 0, the theta-function structure is trivial: 
'd{z,T) oc Y^^5{z + n). For finite disorder, i.e. a > 0, 
the density of the (5-functions is changed making i!){z,T) 
"nearly continuous" a.t T > Ec- The set of m in (fTO|) 
at high temperatures becomes restricted, because only 
the m-configurations directed close to h^^^ contribute to 
the 1?- function. At high temperatures, T > Ec, Eq. 
reduces to 



(11) 



where r 



2tt 



Ta^N 



4S? 



for 



< 1, and X 



eE^^Jn/{2NTEc). If r > 1 then ^{z) « 1 and if r < 1 
then the ?9(2) oscillates with the period 1. While a ^ Ec, 
Ta^N/Ec can be of the order of unity even at weak dis- 
order for T > Ec- Thus, it follows from Eq. ([Til) that the 
disorder-induced term in t dominates and the combined 
action of disorder and temperature makes the current 
a smooth function of voltage. Note that Eq.(fTT|) holds 
provided a > aEc/N\. This condition is not restric- 
tive in large systems and in the thermodynamic limit 
even infinitesimal broadening of quantum levels regular- 
izes I{V), which at high temperatures does not depend 
on the particular choice of a because a chaotic fractalized 
(turbulent) state develops. 

The current-voltage characteristics of the JJA is gov- 
erned by the analytical properties of the i9-function. The 
behaviour of the iJ-function as a function of the dimen- 
sionless voltage and temperature (measured in the units 
of Ec) is illustrated in Fig. [31 One sees that the i?- 
function is smooth at high temperatures but transforms 



into a set of singular peaks at resonant voltages upon 
lowering the temperature. In other words, the 'i?-function 
offers an analytical description of the finite-temperature 
phase transition between conducting and localized states. 



FRACTALIZATION OF THE ENVIRONMENTAL 
SPECTRUM AND UNLOCKED TUNNEL 
TRANSPORT 

The dipoles constituting the environment are com- 
prised of excessive Cooper pairs and Cooper "holes" with 
charges — 2e (+2e) and have the characteristic energy Ec- 
If all the junctions in the JJA were identical, the environ- 
ment energy spectrum would have consisted of narrow 
bands ~ NEj < Ec separated by gaps Ec- Therefore, 
the current would have fiown only at resonant voltages 
V ~ nEc/e, where n is integer. If the charging ener- 
gies are different (disorder in the charging energies also 
accounts for random offset charges that appear in the 
substrate), the tunnelling Cooper pair can relax the en- 



ergy UiEc^'^ to the environment, with rii integer 21 



If the charging energies are incommensurate, the latter 
sum can become arbitrarily small. One can further show 
that a broadening of quantum environmental levels to 
width cr, such that aEc < <7N\, where aEc is the disper- 
sion of the charging energies, turns the local environment 
excitations spectrum effectively quasi-continuous at tem- 
peratures T > Ec- At these temperatures the current- 
voltage (I-V) characteristics is well defined and smooth. 
At T > Ec/ {a^N) the Coulomb blockade modulations in 
the I-V curve almost vanish. If the system is not very 
large such that (for T > Ec) Ta^N < Ec, the modula- 
tions of the I-V curves appear at voltages that are integer 
multiples of Ec, and near T c:^ Ec, the voltage scale of 
current modulations becomes of the order of aEc- 

A quantitative description of fractalization of the envi- 
ronment energy spectrum and the smoothing of Coulomb 
blockade effects can be achieved by expressing the corre- 
lation function as 



IC{t) = (exp[z0(')(i)] exp[-z(/.«(0)]) 
= e ^pjp2...p„ exp 



N 

1=1 



(12) 
(13) 



which is a measure for the temporal evolution 
of the phase determining the current via Is oc 
/ o?t sin(2eyi) Im/C(t). The average is done with re- 
spect to the charging energy part of the Hamiltonian. 
Here ip^"^^ — AE^^h, and the sum is taken over in- 
tegers pi,...,pN, = -2NEc, and Ap^p.^,,,p^ = 



Z-ie-^.2£;(''p?/T^ Noticing that Eq. ^ in fluid dy- 
namics describes the velocity field in the liquid in which 
a turbulent flow develops, one realizes that the time 
evolution of IC(t) is governe d by the Landau-Hopf sce- 
nario of turbulence 0, SEl, 24 1 with Reynolds number 



5 




V/E^ V/E^ V/E^ 

FIG. 3. Current fractalization I-V characteristics at temperatures T/Ec = 10"'^^ [n = 2, . . . , —2, panels a)-e)] according to 
Eq. @. At larger temperatures the /-^/characteristics is approximated by the scalar theta function with a = 0.056 (red dashed 
line). For the highest temperature (n = 2), the slope is given by 2eT//re"^-^'=/^^ (~ 0.13eV/Ec for the used parameters). AU 
plots are for N — S and a — 10~''-Ec. 



TZe = T/Ec- A sequence of incommensurate frequencies 
is generated by bifurcations occurring at TZe = 1. 

The effective level broadening in the high- 
temperature regime, T > i?c, is then 5E = 
2TTEc[a^/i4:EfJ + Ta^N/Ec]- We see that a ther- 
modynamically infinitesimal broadening a of quantum 
levels generates a finite broadening 6E ~ 2'iTTa'^N due 
to the combined action of temperature and disorder 
giving rise to a smooth I-V characteristics. This 
is in complete analogy to introducing a finite width 
of quantum levels of a given subsystem of a larger 
thermodynamic system when proving the equivalence 
of microcanonical, canonical, and grand canonical 
ensembles in quantum statistics. The finite width 
appears due to interaction of this subsystem with the 
thermostat but disappears from the final expressions. 
To conclude at this point, we have demonstrated that at 
high temperatures, T > Ec, the environmental spectrum 
is quasicontinuous and the system is in a conducting 
state. At T < Ec the environmental spectrum retains 
its discrete character and transport is suppressed. The 
localization-dclocalization transition temperature is 
well defined by the condition TZe ~ 1, where the first 
Landau-Hopf bifurcations appears. 

The evolution of the low- voltage I-V curves (|9]) upon 
going from high, T > Ec, to low, T < Ec, temperatures is 
illustrated in Fig. [3l Panels a)-c) show first a practically 
smooth linear I-V dependence at T = lOEc which then 
develops pronounced modulations at resonance voltages 
upon cooling, and finally at T ~ Ec transforms into a 



dense set of resonant spikes. Panels d) and e) show that 
&t T < Ec the Is{V) dependence is a palisade of dis- 
tinct resonant voltages apparently having a hierarchical 
structure, reflecting the hierarchical glassy-like structure 
of the environmental spectrum. 

Having derived the I-V curves at various tempera- 
tures, one can determine the temperature dependence of 
the linear resistance i?(T) at small F — >■ voltages by in- 
tegration of the current over a small voltage interval. At 
high temperatures, T > Ec, the resistance is determined 
by the linear fit of the I-V curve at low voltages. Below 
T = Ec the fractalized hierarchical peak structure of the 
I-V dependence develops and the notion of a smooth I-V 
characteristics ceases to exist. In order to extract a mean- 
ingful resistance value at low temperatures, one finds the 
electric power integrals over a small voltage interval near 
V = of the I — V characteristics and determines the 
resistance from the relation Vj^/R = /g I{V)dV. The 
value Vm is small as compared to Ec/e, but still large 
enough to capture several peaks. Since one can expect 
that the peaks are rather narrow in experiments, the cur- 
rent baseline, i.e. the straight current line limiting the 
current from below, gives a good approximation for the 
resistance as well. 

Asymptotic expressions for the low- and high- 
temperature limits of the resistance can be written as 

i?asym « cxp + St \og{T/Ec)^ , (14) 

with 5t = ~{N - 3)/2 for T/Ec < 1 and St = 3/2 for 
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FIG. 4. Superactivation Dependence of the resistance R on 
inverse temperature using the same parameters as in Fig. [3] 
At temperatures just above Ec Arrhenius behavior is found, 
while at very high temperatures the resistance goes over in a 
power-law. The dotted line separated the linear behaviour of 
the current (delocalized region) and the fractalized behaviour 
(localized state at low temperatures) . Since the I-V character- 
istics becomes fractalized at low temperatures, the resistance 
is extracted by integration of the current at low voltages. 



T/Ec > 1 is derived from the exact Eq. ©. The nu- 
merically computed plot of lni?(r) vs. l/T is shown 
in Fig. m The high temperature limit suggest that 
Arrhenius behaviour should be found for temperatures 
Ec<T < {N/3 - l)EcW~^{N/3 - 1) (for iV > 3, where 
W{x) is the Lambert VF- function, defined as the inverse 
function of f{W) = We^, with W-^{5/3) w 1.3 for 
= 8). Whereas in the low temperature limit, T < E^ 
the deviation from the Arrhenius exponent grows loga- 
rithmically with lowering the temperature. An interme- 
diate Arrhenius behaviour is indicated by a linear fit in 
Fig.m Note, that the 'activation' energy Ea — NEc', this 
reflects the macroscopic character of the Coulomb block- 
ade effect that governs the insulating behaviour of the 
JJA: an extra Cooper pair placed in a granule polarizes 
the whole system, and since the Coulomb energy of two 
interacting charges grows linearly with the separation be- 
tween the charges in ID, the activation Coulomb barrier 
is indeed cx NEr 
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26[. At large temperatures one ob- 
serves an upturn in the lni?(l/r) dependence (see Fig.H]) 
in a nice accordance with the experimental observations 
in TiN films 



DISCUSSION AND CONCLUSIONS 

The derived I-V characteristic for a ID system of 
Coulomb-interacting bosons in the presence of relatively 
strong disorder represents an analytical description of 
a finite temperature phase transition between two dis- 
tinct dynamic states with different conducting and en- 
ergy relaxation properties. In the energy space, this 
transition manifests itself as a transition between a 'lam- 
inar' and 'turbulent' spectral fiow of the environment 




FIG. 5. Resistance at small T Low voltage behavior of 
the current at low temperatures (here T — 0.5Ec) for N — 8, 
a = 10~^Ec, and a — 0.056. This plot shows the hierarchy of 
current peaks (green dotted and dashed lines) and the current 
baseline (red solid line). The current baseline (which limits 
Is/T from below) can be used to extract an approximative 
value for the resistance which is observed in experiments at 
low temperatures. However, a more accurate value is obtained 
by calculating the electric power (see text). 



energy levels occurring at the critical Reynolds number 
7^e* = {T/E,/) = 1. The evolution of the I-V behaviour 
shown in Fig. [3] indicates that for T < Ec the distribution 
of the resonant voltages, i.e. the distribution of charac- 
teristic energies of the environment acquires a hierarchi- 
cal structure (see Fig. [5]). We thus can conjecture that 
at the transition temperature (T* ~ Ec) the dipole envi- 
ronment freezes into a glassy state where dipoles get lo- 
calized, their local spectrum acquires a gap Ec, and they 
become inefficient to mediate energy relaxation of tun- 
nelling carriers. This suggests the one can relate the for- 
mation of a glassy state - which basically can be viewed 
as a transition between the two dynamic states, localized 
and conducting - on a generic level with the onset of 
turbulence in the spectral flow of the respective energy 
spectrum of the system involved. 

The obtained results are valid, strictly speaking, at low 
voltages eV < Ec, since the use of an equilibrium dis- 
tribution function presumes that the time between the 
charge tunnelling events has to be greater than the re- 
laxation time and thus the voltage has to be sufficiently 
low. The observation of a dc tunnelling transport in large 
JJAs described above is feasible in experiments with ar- 
tificially manufactured ID systems with large enough ra- 
tio of the capacitance to the ground of the supercon- 
ducting grains with respect to the average capacitance of 
the junction (C/Co) such that the electrostatic screen- 
ing length A exceeded the perimeter of the system, L. 
Even in the case where A < L one can expect significant 
suppression of the conductivity in the low-temperature 
region, although the characteristic energy scale will re- 
duce to EcN in that case, where N = NA/L, effectively 
replacing the number of junctions by N. Our findings 
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are in accord with early results by 28| where the sup- 
pression of the conductivity at low-temperatures and the 
corresponding voltage threshold behaviour was observed 
in the long ID JJA arrays. Furthermore, our analytical 
results corroborate the findings of [l^ [2^ about the for- 
mation of the superinsulating state with the drastically 
suppressed conductivity a.t T < Ec- 

The results obtained for a JJA ring can be general- 
ized onto 2D JJAs, since the resulting transport across 
2D JJAs can be presented as a superposition of ring- 
like currents over overlapping closed loops spanning the 
2D system. This opens a remarkable opportunity for 
constructing a comprehensive quantitative description of 
the glass formation and localization transition in higher 
dimensions which will be discussed elsewhere. 
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